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ABSTRACT 

ASAS is a long term project to monitor bright variable stars over the whole sky. It 
has discovered 50,122 variables brighter than V < 14 mag south of declination +28°, 
and among them 11,099 eclipsing binaries. We present a preliminary analysis of 5,384 
contact, 2,957 semi-detached, and 2,758 detached systems. The statistics of the distri- 
bution provides a qualitative confirmation of decades old idea of Flannery and Lucy 
that W UMa type binaries evolve through a series of relaxation oscillations: ASAS 
finds comparable number of contact and semidetached systems. 

The most surprising result is a very small number of detached eclipsing binaries 
with periods P < 1 day, the systems believed to be the progenitors of W UMa stars. 
As many (perhaps all) contact binaries have companions, there is a possibility that 
some were formed in a Kozai cycle, as suggested by Eggleton and his associates. 
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1 INTRODUCTION TO ASAS 

ASAS - All Sky Automated Survey, is a long term project 
dedicated to detection and monitoring variability of bright 
stars. This paper presents the results of several years of ob- 
servations done at the Las Campanas Observatory with a 
single instrument: a telescope with the aperture of 7 cm, the 
focal length of 20 cm, done through a standard V-band filter 
and a 2K x 2K CCD camera with 15 fi pixels from Apogee 
(Pojmanski 1997, 1998, 2000, 2002, 2003, Pojmanski and 
Maciejewski 2004, 2005, Pojmanski, Pilecki and Szczygiel 
2005). More information about ASAS is provided on the 
WWW: 

http:/ /www. astrouw.edu.pl/~gp/asas/asas. html 
http : // archive . princeton.edu / " asas/ 
The variable stars were discovered quasi - uniformly for 
declination < +28°, covering almost 3/4 of the full sky. Fig. 
1 shows the distribution of ASAS variables in the sky in the 
Galactic coordinates. The Milky Way is clearly visible, to- 
gether with the dust lanes. The distribution of ASAS stars 
as a function of V-mag is shown on Fig. 2 for all stars 
(17,000,000), all variable stars (50,122), and all eclipsing 
binaries, the latter divided into Eclipsing Contact binaries 
(EC, 5,384), Eclipsing Semi Detached binaries (ESD, 2,957), 
and Eclipsing Detached binaries (ED, 2758). The statistics 
for stars with 8 < V < 12 mag appears to be approximately 
complete, but the efficiency falls rapidly within the range of 
12 < V < 14 mag as the detection limit is approached. Also, 
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the statistics deteriorates for 8 < V because of saturation ef- 
fects. 

All stars were observed for about 5 years, with a small 
subset for 8 years. Typical number of photometric measure- 
ments was several hundred. The distribution of this number 
is shown in Fig. 3. The total number of photometric mea- 
surements was 2,916,000. As ASAS continues its operation 
the number of measurements will increase, approximately 
100 V-band photometric measurements per year per vari- 
able. We intend to continue the project indefinitely, with 
some upgrades. While only V-band results are reported here, 
the I-band photometry was accumulated, and several years 
of data are already stored on a RAID-5 disk system. How- 
ever, it will take another year to process I-band data, well 
over a Tera-byte. There are more stars detectable in the I- 
band, so we expect that the number of variables will more 
than double. We are also planning an expansion of ASAS to 
the northern hemisphere, to fully cover the whole sky. 

The results presented in this paper are not final in any 
sense. Please note, that all data are public domain. A differ- 
ent classification of binaries can be readily done by whoever 
feels like verifying and/or correcting our presentation. 

This is an observational paper, but section 2 gives a 
short introduction to the main ideas about the structure 
and evolution of contact binaries. Section 3 gives the in- 
formation about ASAS data, the classification scheme, and 
some examples of a diversity of light curves. Section 4 pro- 
vides simple statistics of ASAS binaries. Finally, in section 
5 we make a somewhat speculative discussion based on this 
observational paper. 
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Figure 1. The distribution of 50,122 ASAS variables in the sky in Galactic coordinates. The Milky Way is clearly seen, as well as the 
patches of interstellar extinction. The equator is shown with a dashed line. The distribution is limited to declination < +28°. 
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Figure 2. The distribution of ASAS stars as a function of V- 
band magnitude. The total number of stars is about 17,000,000. 
The total number of variable stars is 50,122. The total number 
of eclipsing binaries is 5,384 for contact systems (EC), 2,957 for 
semi-detached systems (ESD), and 2,758 for detached systems 
(ED). Notice, that the efficiency of discovering variable stars de- 
clines for V > 12 because the detection limit is approached, and 
for V < 8 because of saturation effects. 



2 INTRODUCTION TO CONTACT BINARIES 

Contact binaries, also known as W UMa stars, are in a phys- 
ical contact, with continuously changing brightness because 
of large tidal distortion of the two components. 

The first theoretical milestone in the understanding of 
contact binaries was due to Lucy (1968a, b), who proposed 
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Figure 3. The histogram of the number of photometric mea- 
surements obtained during 5 years of ASAS life. A small subset 
of data extends back to 8 years. 



that the two components share common envelope with the 
same entropy, thereby making the effective temperature al- 
most constant over the surface of the two stars. As con- 
tact binaries have the mass ratio distinctly different from 
one, most nuclear energy is generated in the more massive 
component and it is redistributed around the whole surface 
through a moderately thick convective envelope. 

The second theoretical milestone was the recognition 
of the consequences of the fact that the mass - radius rela- 
tion for Zero Age Main Sequence (ZAMS) stars is much 
steeper than for the two Roche lobes. There can be no 
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stable equilibrium between the two stars with a common 
envelope. The system evolves through a sequence of relax- 
ation oscillations, with the mass flowing from star A to B, 
next from B to A, etc. (Flannery 1976, Lucy 1976, Robert- 
son and Eggleton 1977). The cycle repeats on a thermal 
(Kelvin-Helmholtz) time scale. According to thermal relax- 
ation model the binary oscillates between thermal contact, 
with the two eclipses of almost equal depth, and a semide- 
tached phase in which one eclipse is much deeper than an- 
other. 

Hazlehurst (1970) suggested that nuclear evolution of 
the primary component of a contact binary affects its struc- 
ture. Stepieh (2003, 2005) suggested that the currently more 
massive primary was originally the less massive of the two. 
The nuclear burning formed a small helium core, the star 
expanded and transferred mass to the original secondary. In 
analogy with Algol systems the currently more massive com- 
ponent is the less evolved, while the present secondary has 
a small helium core, and it is more advanced in its nuclear 
evolution. 

It is interesting that there is a controversy about ther- 
mal relaxation oscillations in W UMa systems. Some authors 
claim there are no such oscillations (Webbink 2003), while 
others claim that such oscillations exist (Quin 2003, Yakut 
and Eggleton 2005, and references therein). ASAS statistics 
resolves this controversy on purely observational grounds. 



3 ASAS DATA 

Close binaries with a deep common envelope are in thermal 
contact and they have eclipses of almost equal depth. If the 
contact is shallow, or if there is no physical contact, then 
the effective temperatures of the two stars are different, and 
the two eclipses have different depth. 

Theoretical models of relaxation oscillations indicate 
that the radii of the two components change relatively little 
throughout the cycle (Flannery 1976, Lucy 1976, Robertson 
and Eggleton 1977, Yakut and Eggleton 2005, and refer- 
ences therein). With the geometry of the two stars almost 
unchanged, tidal distortions due to geometry remain almost 
the same, and the most profound difference in the light curve 
is the relative depth of the two eclipses. 

For the purpose of this paper the classification of eclips- 
ing binaries was done by decomposing their light curves 
into Fourier coefficients: ai,a2,a3,a4, following Rucihski 
(1997a,b, 1998, 2002) and Pojmahski (2002). The classifi- 
cation is shown in Fig. 4 in the a2 — a4 plane, with the EC, 
ESD, and ED binaries shown with different symbols: ESD 
stars have smaller marks than EC and ED. Some examples 
of bright binaries of EC, ESD, and ED type are shown in 
Fig. 5. The parameters of these 15 binaries are listed in Ta- 
ble 1. The type of an eclipsing binary, EC, ESD, or ED, is 
indicated at the upper right corner of the light curve. The 
complete version of Table 1, for all eclipsing binaries used in 
this paper, and with all eight Fourier parameters, is in the 
file Fourier. E on the web page: 

http://www.astrouw.edu.pl/~gp/asas/asas_paper_data.html 

http://archive.princeton.edu/~asas/asas_paper_data.html 
where all the data on which this paper is based are available 
electronically. These include all photomertic measurements 




a2 

Figure 4. Classification of eclipsing binaries in the Fourier coef- 
ficients plane a2 — a4. The three symbols refer to contact (EC), 
semidetached (ESD), and detached (ED) binaries, following Po- 
jmanski (2002). Note: when the amplitude of variability is very 
small the classification is very uncertain, as in the upper right 
hand corner. 



Table 1. The most important Fourier parameters for light curves 
from Fig. 5 are listed along with orbital periods. Each row is sep- 
arated with a horizontal line. The a,2 parameter can be easily 
translated to the amplitude, a\ and G13 determine the difference 
between the two eclipses, and b\ is related to the difference be- 
tween the two maxima. Full Table 1 is given electronically as file 
Fourier. E on our web page. 



ASAS ID 


P [days] 


H 


"2 


"3 


~4 


bi 


125340-5010.6 
111621-5427.2 
220425-0603.4 


0.4047 
0.4971 
0.7273 


-0.003 
-0.026 
-0.062 


-0.158 
-0.123 
-0.139 


-0.003 
-0.012 
-0.034 


-0.028 
-0.020 
-0.027 


-0.003 
-0.000 
0.021 


181437-1324.0 
125028-6039.8 
162959-4515.5 


1.3992 
2.0439 
1.1376 


-0.001 
-0.030 
-0.068 


-0.110 
-0.085 
-0.155 


-0.004 
-0.010 
-0.032 


-0.054 
-0.031 
-0.049 


0.004 
-0.001 
0.001 


072345-6500.6 
185503-3359.5 
175924-0954.7 


3.1240 
1.1678 
2.3438 


-0.006 
-0.019 
-0.115 


-0.066 
-0.070 
-0.119 


-0.003 
-0.015 
-0.090 


-0.060 
-0.060 
-0.088 


-0.001 
0.002 
0.006 


123140-6237.3 
172116-3759.2 
012917-7243.4 


1.9307 
5.2056 
180.60 


-0.009 
-0.021 
-0.050 


-0.184 
-0.135 
-0.190 


-0.004 
-0.005 
0.016 


-0.051 
-0.033 
-0.052 


0.001 
0.007 
0.041 


051700-5555.4 
105554-4349.6 
132831-4107.9 


0.7904 
1.1541 
16.515 


-0.059 
-0.034 
-0.010 


-0.126 
-0.185 
-0.069 


-0.028 
-0.019 
-0.011 


-0.021 
-0.030 
0.003 


0.020 
0.037 
-0.015 



in compressed files: ec.tgz (49MB), esd.tgz (27MB), ed.tgz 
(27MB). 

Some investigators (e.g. P. Eggleton) strongly pre- 
fer classification of close binaries into EW, EB, and EA 
types, on the grounds that photometry alone cannot pro- 
vide unique geometry. Our classification into EC, ESD, and 
ED types should be considered preliminary, as we have only 
single band photometry and no spectroscopic information 
for thousands of our binaries. 
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Figure 5. Examples of light curves of 9 contact (EC), 3 semi-detached (ESD) and 3 detached binaries (ED) are shown. All these are 
ASAS discoveries, even though all of them are very bright. Notice a distinct difference in the light curve shapes between contact and a 
semi-detached binaries. The first three rows are sorted by minimum depth difference, the fourth row presents contact binaries with long 
end exceptionally long periods, and the fifth with the distinct maximum height difference. See Table 1 for details. 



A sample of the catalog of ASAS eclipsing binaries is 
shown in Table 2. The full table has the information for all 
11,099 binaries, and can be found on our web page, in the 
file Vars.E. 

A subset of all ASAS results related to binary stars was 
'frozen' for the epoch of this paper. ASAS is an on-going 
project, with more data added every clear night. Therefore, 
as time goes on, the volume of data related to eclipsing bi- 
naries will increase, and the quality will improve. To make 
the results presented in this paper reproducible we decided 
to provide 'frozen' data. 

Note: all light curves presented in Fig. 5 are for bright 
variables discovered with ASAS survey. They were not 
known before. Just as expected, the majority of ASAS vari- 
ables are new discoveries (Paczynski 1997). 

4 STATISTICS OF CONTACT BINARIES 

All the statistics of this section is based on data provided in 
the file Vars.E, which can be found on our web page. Table 



2 is a sample of this file. In particular, the classification of 
ASAS binaries is given there. 

The Galactic distribution of ASAS eclipsing binaries 
shows a significant difference, with detached binaries most 
strongly concentrated to the galactic plane, while the short 
period contact binaries have almost isotropic distribution. 
It is well known that there is a period - luminosity relation 
for W UMa stars (Rucinski 1996, Klagyivik and Csizmadia 
2004), and the Galactic distribution of those stars is not 
surprising. The strong concentration of detached binaries to 
the Galactic plane implies these are even more luminous and 
massive, on average. 

The binary period distribution is shown in Fig. 6 for 
stars with the Galactic latitude |b| > 30° to reduce the lumi- 
nosity bias. The distribution of contact systems (EC) peaks 
near 0.37 days, it has a sharp cut-off at 0.2 days, and a long 
tail extending beyond 1 day. In fact the tail extends to more 
than one hundred days, as it is apparent in Fig. 8 and Fig. 
10. Also shown is the distribution of orbital periods of semi 
- detached (ESD) and detached (ED) binaries. 

The important topic ot this paper is the presentation 



Eclipsing binaries in ASAS catalog 5 



Table 2. A sample of the catalog of ASAS eclipsing binaries is presented here. Each star has its own ASAS designation 
(ID), equatorial coordinates (RA and DEC), orbital period, epoch of minimum light (TO), and ASAS classification. 
Additional information of other designation and classification is also given, when available. Full Table 2 is available 
electronically as the file Vars . E in our web page. 

ID RA DEC Period TO V Amp Class Other Other 

(2000) (2000) [days] 2450000+ [mag] [mag] - ID Class 



065819+1028.4 


06:58:19 


+ 10:28:24 


0.447763 


2387.63 


11.04 


0.28 


EC 






070017+0202.2 


07:00:17 


+02:02:12 


1.397550 


1873.07 


10.53 


0.57 


EC 


V0460 


Mon EB/KE 


211201-2003.2 


21:12:01 


-20:03:12 


0.352202 


1872.84 


13.22 


0.34 


EC/ESD 






025535-0219.9 


02:55:35 


-02:19:54 


0.79274 


1920.00 


11.59 


0.54 


ESD 






095930-6440.1 


09:59:30 


-64:40:06 


2.2428 


1869.42 


11.48 


1.12 


ESD 






202110-4333.9 


20:21:10 


-43:33:54 


1.00218 


1876.053 


12.08 


0.71 


ESD/ED 


V2265 


Sgr EA 


001855-7954.9 


00:18:55 


-79:54:54 


0.90310 


1870.45 


11.28 


0.77 


ED 






114500-7247.6 


11:45:00 


-72:47:36 


22.23037 


1895.60 


11.15 


0.09 


ED 






235052-2316.7 


23:50:52 


-23:16:42 


1.4023 


1871.06 


9.42 


0.29 


ED 
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Figure 6. The distribution of periods of ASAS contact binaries 
(EC) at high Galactic latitude |b| > 30° is plotted. The distribu- 
tion peaks near 0.37 days, it has a sharp cut-off at 0.2 days, and 
a tail extending far beyond 1 day. Also shown is the distribution 
of orbital periods of semi - detached (ESD) and detached (ED) 
binaries. Contact binaries outnumber other binaries for binary 
periods shorter than 1 day. 



of eclipse depth of EC binaries, as this is directly related to 
the issue of relaxation oscillations. For contact binaries the 
depth of both eclipses is directly related to the Fourier coef- 
ficients of the light curve. The following analytical formulae 
approximate these relations with better then 5% accuracy 
for the majority of stars: 



(l.2a 2 - 2a2 + (aj + a 3 )) 



D s = -2 (l.2a 2 - 2a^ - (aj + as)) 



[magj . 



[magj 



(1) 
(2) 



correspond to the depth of the primary 
and secondary minimum, respectively. Three examples are 
presented in Fig. 7. 

The distribution of primary eclipse depth with binary 
period is shown in Fig. 8. A histogram of primary eclipse 
depth integrated over binary period is shown in Fig. 9. Note 
that while the vast majority of W UMa binaries have periods 
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Figure 7. The first row of EC stars from Fig. 5 is shown. Ver- 
tical lines represent depth of a primary (D p ) and a secondary 
(D s ) eclipse calculated with equations (1) and (2), respectively. 
Flat bottom eclipses may result in a significant error, which is 
somewhat reduced in the D s /D p ratio. 



in the range 0.2 < P < 1.2 days, there are contact binaries 
with periods over 100 days, as it is apparent in Fig. 5 and 
Fig. 8. 

The eclipse depth ratio is shown in Fig. 10 for EC bi- 
naries, and in Fig. 11 for EC and ESD binaries. These are 
important figures. They demonstrate that some EC bina- 
ries are in good thermal contact, as the depth of secondary 
eclipse is almost the same as the primary. But many other 
EC binaries have the two eclipses which are very unequal, 
implying that they are not in thermal contact, yet their ge- 
ometry is close to that expected for a contact system. Fi- 
nally, there are ESD binaries, in which the geometry is very 
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Figure 8. The distribution of primary eclipse depth as a function 
of orbital period for contact binaries. D p is the fraction of light 
obscured in the primary eclipse. 



Figure 10. The ratio of eclipse depth is shown for contact bina- 
ries as a function of their period. Stars in thermal contact have 
the eclipse depth ratio close to one. 
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Figure 9. The distribution of primary eclipse depth integrated 
over all periods. 



different, and of course the two eclipses have usually differ- 
ent depth. Qualitatively, this is just what was expected in 
the models with relaxation oscillations (cf. Lucy 1976, Flan- 
nery 1976, Robertson and Eggleton 1977, Yakut and Eggle- 
ton 2005). There is no shortage of stars without thermal 
contact. Combining ASAS data with models of relaxation 
oscillations should provide quantitative verification of the 
theory, but this task is beyond the scope of our paper. 

Contact binaries often have maxima of different height, 
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Figure 11. The ratio of the eclipse depths, integrated over bi- 
nary periods, is shown for contact (EC) and semi-detached (ESD) 
systems. Stars in a thermal contact have the eclipse depth ratio 
close to one. Note the step in the histogram at D s /D p « 0.7 The 
same feature is also apparent in Fig. 10 at logP f» —0.4. 



with the maximum following the primary (i.e. deeper) 
eclipse being either higher or lower than the maximum fol- 
lowing the secondary eclipse. The distribution of the differ- 
ence is shown in Fig. 12. The positive value of the bl/|a2| pa- 
rameter indicates that the maximum following the primary 
eclipse is brighter of the two. The asymmetry in distribution, 
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Figure 12. The distribution of light curve asymmetry of contact 
binaries is shown - this is known as O'Connell (1951) effect. Most 
stars have their light maxima of approximately equal height, but 
there is an asymmetry in the distribution: the maxima that follow 
the primary minimum are on average higher than the maxima 
preceding the primary eclipse. The O'Connell effect is likely a 
consequence of gas streams in these binaries. 



i.e. the excess of positive values of bl/|a2|, is known as the 
O'Connell (1951) effect (cf. also Rucihski 1997b): the max- 
imum following the primary eclipse is on average brighter, 
presumably due to stream of gas flowing between the two 
stars. An extreme case of this phenomenon is shown in V361 
Lyr (Kaluzny 1991). 



5 DISCUSSION 

Our conclusion, based on the distribution of eclipse depths 
(Fig. 10, 11), is that relaxation oscillations, first proposed by 
Lucy and Flannery, are real. There is no shortage of binaries 
corresponding to no thermal contact, with a very different 
depth of the two eclipses. Model calculations of the type 
recently done by Yakut and Eggleton (2005) , when combined 
with ASAS data, should allow a quantitative verification of 
the theory. 

A very large number of contact or near contact sys- 
tems reveal features never noticed before: the distribution 
of eclipse depth ratios as shown in Fig. 10 and Fig. 11 has 
a distinct break around logP « —0.4 and D s /D p w 0.7. We 
do not speculate on the origin of this feature, but we bring 
this break to the attention of our readers. This break is best 
seen in systems with the most robust classification. 

The most surprising result of this paper is presented in 
Fig. 6. A traditional view for the origin of W UMa contact 
binaries is to assume that they come from detached binaries 
of comparable periods. For the first time we have approxi- 
mately complete statistics of binaries of all types with or- 
bital periods shorter than one day, and there are very few 
detached binaries. Obviously, they are more difficult to find 



than either contact or semi - detached systems. As times 
goes on the statistics of ASAS detached binaries will im- 
prove, and a statistical analysis will tell us if there is a prob- 
lem with the origin of W UMa stars. At this time the contact 
systems seem to appear "out of nowhere". 

In about a year we shall have I-band data for our eclips- 
ing binaries. This will make it possible to be more quantita- 
tive about the distribution properties, including the distri- 
bution in our Galaxy. Also, it will be much easier to quantify 
the impression about a shortage of detached binaries and the 
space density of systems of different types: EC, ESD, and 
ED. At this time it is premature for us to speculate about 
the outcome of binary statistics while we wait for the I-band 
data. 

While this is an observational paper, written to promote 
the usefulness of the ASAS catalog of variable stars, we are 
tempted to speculate about possible interpretation of Fig. 
6, which was so surprising to us. We are not consistent with 
the previous paragraph, but the temptation is hard to resist. 
The following is a speculative hypothesis. 

There has been a gradual emergence of the notion that 
contact systems have companions (Rucihski and Kaluzny 
1982, Chambliss 1992, Hendry and Mochnacki 1998, also 
Tokovinin 2004). Recently. Pribulla and Rucihski (2005) 
found that up to 50% of W UMa binaries have companions. 
This opens up a possibility that Kozai (1962) cycle operates 
in some such triples, as suggested by Kiseleva, Eggleton, 
and Mikkola (1998), and by Eggleton and Kiseleva-Eggleton 
(2001). The mechanical three-body orbital evolution, with a 
large range of inclinations of the two orbits, inner and outer, 
may induce large variations in the eccentricity of the inner 
binary. In time the eccentricity may become large enough to 
make a contact system out of the inner binary. 

In this scenario the inner binary has the initial period 
that is relatively long, and in most cases no eclipses would 
be detectable. During a Kozai cycle the inner binary occa- 
sionally has a chance to reach a physical contact: it may 
either merge forming a single star, or it may become a W 
UMa star. This process is somewhat similar to the model of 
formation of close binaries in globular clusters (Pooley et al. 
2003, and references therein). 

The surprisingly small number of short period detached 
eclipsing binaries as seen in Fig. 6 may indicate that the 
Kozai cycle is not just a curiosity, but it may be an impor- 
tant channel for forming W UMa stars. This possibility can- 
not be rejected without careful analysis. After all W UMa 
stars are rare, with the local space density of just 0.2% of the 
main sequence stars (Rucinski 2002, and references therein). 



ACKNOWLEDGEMENTS 

We are very grateful to P. P. Eggleton, D. Fabrycky, S. 
Rucihski and K. Stepieri for many helpful discussions. 

The ASAS project was supported to a large extent by 
the generous donation of William Golden. This work was 
partially supported by the Polish MNil grant no. 1P03D- 
008-27. 



REFERENCES 

Chambliss, C. R. 1992, PASP, 104, 663 



8 B. Paczynski, D. Szczygiel, B. Pilecki, and G. Pojmanski 



Eggleton, P. P., Kiselcva, L. 2001, ApJ, 562, 1012 

Flannery, B. P. 1976, ApJ, 205, 217 

Hazlehurst, J. 1970, MNRAS, 149, 129 

Hendry, P. D., Mochnacki, S. W. 1998, ApJ, 504, 978 

Kaluzny, J. 1991, AcA, 41, 17 

Klagyivik, P., and Csizmadia, Sz. 2004, PADEU, 14, 303 

Kiseleva, L. G., Eggleton, P. P., Mikkola, 1998, MNRAS, 300, 292 

Kozai, Y. 1962, AJ, 67, 579 

Lucy, L. B. 1968a, ApJ, 151, 1123 

Lucy, L. B. 1968b, ApJ, 153, 877 

Lucy, L. B. 1976, ApJ, 205, 208 

O'Connell, D. J. K. 1951, Riverview Pub. 2, 85 

Paczynski, B. 1997, vsar.conf, 357 

Pojmanski, G. 1997, AcA, 47, 467 

Pojmanski, G. 1998, AcA, 48, 35 

Pojmanski, G. 2000, AcA, 50, 177 

Pojmanski, G. 2002, AcA, 52, 397 

Pojmanski, G. 2003, AcA, 53, 341 

Pojmanski, G., and Maciejewski, G. 2004, AcA, 54, 153 
Pojmanski, G., and Maciejewski, G. 2005, AcA, 55, 97 
Pojmanski, G., Pilecki, B., Szczygiel, D. 2005, astro-ph/050817 
Pooley, D. et al. 2003, ApJL, 591, L131 
Pribulla, T., Rucinski, S. M. 2005, private communication 
Quin, S. 2003, MNRAS, 342, 1260 

Robertson, J. A., and Eggleton, P. P. 1977, MNRAS, 179, 359 

Rucinski, S. M. 1996, ASPC, 90, 270 

Rucinski, S. M. 1997a, AJ, 113, 407 

Rucinski, S. M. 1997b, AJ, 113, 1112 

Rucinski, S. M. 1998, AJ, 115, 1135 

Rucinski, S. 2002, PASP, 114, 1124 

Rucinski, S. M., Kaluzny, J. 1982, Ap&SS, 88, 433 

Stepien, K. 2003, IAUS 219, 134 

Stepien, K. 2005, astro-ph/0510464, submitted to A&A 

Tokovinin, A. 2004, RMxAC, 21, 7 

Webbink, R. F. 2003, ASPC, 293, 76 

Yakut, K., Eggleton, P. P. 2005, ApJ, 629, 1055 



This paper has been typeset from a TffjX/ WFftK file prepared 
by the author. 



